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Abstract. We present atmospheric parameters and Fe abundances de- 
rived for the majority of dwarf stars (north of -30 degrees declination) 
which are up to now known to host extrasolar planets. High-resolution 
spectra have been obtained with the Sandiford Echelle spectrograph on 
the 2.1m telescope at the University of Texas McDonald Observatory. 
We have used the same model atmospheres, atomic data and equivalent 
width modeling program for the analysis of all stars. Abundances have 
been derived differentially to the Sun, using a solar spectrum obtained 
with Callisto as the reflector with the same instrumentation. A similar 
analysis has been performed for a sample of stars for which radial velocity 
data exclude the presence of a close-in giant planetary companion. The 
results are compared to the recent studies found in the literature. 



1. Introduction 

To examine the relation between the metallicity of the stellar host and the 
existence of close-in massive planets, we have analyzed high-resolution, high 
signal-to-noise spectra of the dwarf stars with super-Jupiter planets listed in 
the "Extrasolar Planets Encyclopaedia" 1 (44 stars at the time when this poster 
was prepared, hereafter CGP dwarfs). Note that this list does not include stars 
with companions with M p sini > 12Mj. A list of comparison stars was compiled 
using the results of the Lick planet search (Cumming et al. 1999). We analyzed 
all stars for which radial velocity data exclude the presence of a close-in giant 
planetary companion (23 stars, hereafter no-CGP dwarfs). The selection process 
is illustrated in Figure 1, which shows the upper limits for planetary masses (in 
Jupiter-masses) as a function of orbital radii (in AU) for this sample, as derived 
from radial- velocity data, and masses and semi-major axes for known extrasolar 
planets. Those of the latter which would fall below the lower limit for the 
comparison sample have been excluded. 

2. Observations 

The spectra were obtained at the 2.1m telescope of the McDonald Observatory 
with the Sandiford Echelle spectrograph. They cover a wavelength range from 



http://www.obspm.fr/encycl/encycl.html 
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Figure 1. Red dots: Masses and semi-major axes of currently known 
extrasolar planets. Blue lines: Upper (mass) and lower (orbital radius) 
limits for stars from the Lick survey (Cumming et al.1999). 
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about 4840 to 7000 A with a resolution of 60000. The list of the observed stars 
can be found in Tables 1 and 2. The reductions were done with IRAF (echelle 
order extraction) and a Windows based graphical package (ASP) for continuum 
and wavelength setting and equivalent width determination developed by REL. 
A comparison of equivalent widths with published values indicates agreement at 
the 5% level. 

3. Analysis 

The abundance analysis was done strictly differentially with respect to the Sun 
using a large sample of lines (cf. Figures 4a, b). We have obtained a solar flux 
spectrum (Callisto) using the same instrumentation and reduction procedure as 
used for our program stars. Figure 2 shows a portion of the observed spectra of 
HD 195019 and the Sun, as well as a synthetic spectrum, computed with SYNTH 
(Piskunov 1992). Synthetic equivalent widths were calculated and fit to the 
observations by variation of abundance for each line with the program LINES 
(originally by Sneden 1974). We used MARCS model atmospheres (Gustafsson 
et al. 1975) and an atomic data list compiled by REL. The relative abundance 
for each line was then calculated by subtracting the abundance obtained from 
the Callisto spectrum for the corresponding line. From all Fe lines available 
within the observed spectral range, only those giving an abundance within two 
standard deviations (2a) from the mean have been retained for the derivation 
of the atmospheric parameters described in the following subsection. This was 
done in order to discard "outliers" automatically. For the calculation of the final 
mean abundances, however, the line list has been further reduced by rejecting 
all lines with an abundance deviation from the mean of more than la, which 
decreased the errors of the abundances significantly. 

3.1. Atmospheric parameters 

To determine the atmospheric parameters, iron line abundances were calculated 
for each star for a small (T e g, logg, v m i cr o) ~~ grid, centered on t> m icro = 1 kms -1 
and (T e g, logg) determined by one of three methods in the following order of 
preference: 

1. Literature data (previous abundance analyses 2 ) 

2. Calibration for Geneva photometry (Kiinzli et al. 1997) 

3. A new calibration for line-depth ratios and T e g for our sample of dwarf 
stars. We used 20 of the 32 line-pairs which Kovtyukh and Gorlova (2000) 
had used for a similar calibration for supergiants, our observations, and 
effective temperatures from methods 1 and 2. An example for a useful 
(left panel) and a rejected (right panel) line-pair is shown in Figure 3. 
The line-ratio - T e g- relations were obtained from polynomial fits to the 
data, or calculated from synthetic spectra in some cases where only few 
datapoints were available. 



2 Santos et al. (2000, 2001), Gonzalez et al. (1997-2001), Fuhrmann (1998), Edvardsson et al. 
(1993) 
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Table 1. Observations of CGP dwarfs. 



rilJ 


date 


XlJJ 


date 


riJJ 


date 


8574 


2001-08-23 




1997-02-03 


169830 


2000-08-18 




2001-08-27 




1999-10-19 




2000-08-21 


9826 


1999-10-19 


80606 


2002-01-23 


177830 


2000-08-18 




1999-10-23 




2002-01-28 




2000-08-21 


10697 


1992-11-03 


82943 


2002-01-22 


178911 


2001-05-08 




1992-11-05 




2002-01-25 




2001-05-14 




1992-11-07 


92788 


2001-05-08 




2001-08-27 




1997-02-03 




2001-05-14 


179949 


2001-08-24 




1999-10-19 


89744 


2001-05-08 




2001-08-28 


12661 


2000-08-18 




2001-05-14 


186427 


1999-05-26 




2000-08-21 


106252 


2001-05-08 




1999-05-28 


19994 


2000-10-20 




2001-05-14 




2000-08-19 




2001-02-11 


114762 


1999-05-26 


187123 


1999-05-30 




2001-08-23 




1999-05-28 




1999-10-25 




2001-08-25 




2000-01-25 


190228 


2000-10-22 


28185 


2001-10-25 


117176 


1999-05-25 




2001-08-23 




2001-10-27 




1999-05-28 




2001-08-27 


33636 


2002-01-23 




2000-01-25 


192263 


2000-08-19 




2002-01-26 


120136 


1999-05-25 




2000-08-21 


37124 


2000-01-26 




1999-05-28 


195019 


1999-05-30 




2000-01-29 




2000-01-25 




1999-10-25 


38529 


1992-11-04 


130322 


2000-01-26 


209458 


2000-08-18 




1992-11-05 




2000-01-30 




2000-08-21 




1992-11-07 


134987 


2000-04-25 


210277 


1999-10-19 




1993-03-10 




2000-04-30 




1999-10-25 




1997-01-31 


136118 


2002-05-21 


217014 


1992-11-03 




1997-02-02 




2002-05-24 




1992-11-06 


50554 


2001-10-23 


141937 


2001-05-08 




1992-11-07 




2001-10-27 




2002-01-26 




1997-10-18 


52265 


2001-10-24 


143761 


1999-05-26 




1999-10-19 




2001-10-27 




1999-05-28 


217107 


1999-10-19 


68988 


2002-01-24 




2001-05-10 




1999-10-24 




2002-01-26 


145675 


1993-03-06 


222582 


2000-08-18 


74156 


2001-10-25 




1993-03-09 




2000-08-21 




2001-10-28 




1993-03-11 


BD -10 3166 


2001-05-10 


75732 


1993-03-03 




1997-04-24 




2002-01-26 




1993-03-05 




1998-05-15 




2002-01-28 




1993-03-11 


168443 


1999-05-30 


Callisto 


2001-10-26 




1993-03-07 




1999-10-25 




2001-10-28 
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Table 2. Observations of no-CGP dwarfs. 



HD 


date 


HD 


date 


HD 


date 




1 999-1 (1-9^ 


39147 


1 999-1 1-03 


14Qfifi1 


1 993-03-03 




9nnn n& 99 




1 QQ9 1 1 0^ 






4628 


1999-10-25 




1992-11-08 




1993-03-07 




2000-08-20 




1993-03-10 




1993-03-10 




2000-08-22 




1995-10-11 




1997-04-24 


10476 


2000-08-20 




1997-01-31 




1997-04-28 




2001-08-23 




1997-02-02 


157214 


2000-04-25 


12235 


1992-11-04 


48682 


2000-10-20 




2000-04-29 




1992-11-06 




2001-10-25 


170657 


2001-08-24 




1992-11-08 




2001-10-27 




2001-08-27 




1992-11-09 


50281 


2001-10-26 


166620 


2000-10-19 




1997-02-03 




2002-01-25 




2001-08-25 




1999-10-19 


76151 


1993-03-04 




2001-08-27 


16160 


2000-08-18 




1993-03-06 


185144 


1999-10-19 




2001-08-25 




1993-03-07 




1999-10-23 


16895 


2000-08-19 




1993-03-10 


186408 


1999-05-26 




2001-10-23 




1997-02-03 




1999-05-28 


22484 


1999-10-21 




2000-01-29 




2000-08-19 




1999-10-24 


84737 


2002-05-21 


201091 


1997-08-17 


26965 


1999-10-19 




2002-05-24 




1999-10-19 




1999-10-23 


126053 


2002-05-21 


219134 


1999-10-19 








2002-05-24 




1999-10-24 










222368 


2000-08-19 



2000-08-21 
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Figure 2. Sample spectra of HD 195019 (green dots) and the Sun 
(Callisto, blue dots). Red line: synthetic spectrum. Note that for the 
Fe I line at 5317 A the gf value is obviously wrong, but this source of 
error is canceled out in a strictly differential analysis. 
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Figure 3. Two examples for the dependence of line-depth ratios on 
effective temperature. Red: literature data, blue: Geneva photometry, 
green: polynomial fit, gray: relation for ratios calculated from synthetic 
spectra. The relation shown in the right panel has not been used for 
the calibration. 



The final atmospheric parameters were determined by driving the slopes of 
the line-strength and excitation energy relations with Fe I line abundance to zero 
and demanding ionization equilibrium. As an example, Figure 4 shows the trends 
visible in an abundance versus equivalent width or excitation energy diagram, 
if T e Q, logg, and f m icro deviate from the best-fit values by 100 K, 0.1 [cgs], 
and 0.3 kms -1 , respectively (for the planet host HR 5072). We determine the 
best-fit values by first identifying the value of v m i cro , for which the trends in 
the above mentioned diagrams are zero, for each model in the (T e s, log g) grid 
(cf. Figure 5). Figure 6 shows these "zero-points" for equivalent width (dots) 
and excitation energy (diamonds) for each grid point. From the set of (T e g, 
log g) models for which the two zero-points are equal (i.e. the intersections 
of two corresponding lines in Figure 6), we can finally identify one parameter 
pair, which gives the same abundance for neutral and ionized Fe lines. Figure 7 
illustrates the dependence of the abundance difference [Fe I] — [Fe II] on T e g and 
log g. The dependence of this quantity on v m \ cro is negligible. The combination 




Figure 4. Fe abundances in dex (0.00. . . [Fe]Q) versus atomic data for 
HR 5072 and (T e s, logg, w m icro) = (5700, 4.0, 0.7). These parameters 
obviously do not correspond to the best-fit values. 
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Figure 5. Dependence of the slopes appearing in diagrams like those 
shown in Figure 4 on v m i cro and T e g, for Fe I, here for logg=4.1 (HR 
5072). The errorbars are determined from linear regression. 
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Figure 6. Values of i> m icro) for which the slopes in the abundance ver- 
sus atomic data diagrams vanish, in dependence of T e g- and log g. Dots 
and diamonds correspond to equivalent width and excitation energy, 
respectively. 




T eff [K] 
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5500 



4.0 4.1 

log g [log cm/s 2 ] 

• 5600 



4.4 



5700 
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Figure 7. Difference between the mean abundances of neutral and 
ionized iron, for v m i cro =1.0 km s -1 . The errorbars correspond to the 
standard deviation of the line abundances. 
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Figure 8. Red: Intersections of two lines from Figure 6 corresponding 
to equivalent width and excitation energy. Blue: Combinations of T e g 
and log<7, for which the difference between the mean abundances of 
neutral and ionized iron is zero (determined from Figure 7). 
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Figure 9. Variation of the mean iron abundance with parameters. 
Here, the error bars indicate the abundance range covered by using the 
iWro range of Figure 5. 
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of information extracted from Figures 6 and 7 is shown in Figure 8. Figure 9 
shows the variation of the mean Fe abundance within the parameter grid. From 
these diagrams we estimate the errors in our final spectroscopically determined 
parameters to be 50 K, 0.05 [cgs] and 0.05 kms" 1 , for T c g, logg, and v m i CIO , 
respectively. 

4. Results 

4.1. Atmospheric parameters 

A comparison of the parameters derived as described above with parameters 
found in the literature and derived from photometry results in the mean differ- 
ences and standard deviations listed in Table 3 (parameter [this work] — param- 
eter [comparison] ) . 

Table 3. Differences between parameters determined in this work, in 
the literature and with photometric calibrations. 





N 


T e fi 




log 


9 


[Fe] 




Spectroscopy 1 


52* 


55 


± 


90 


0.05 


± 


0.20 


0.02 


± 


0.07 


CGP dwarfs 


37 


50 


± 


90 


0.05 


± 


0.20 


0.02 


± 


0.08 


No-CGP dwarfs 


9 


80 


±120 


0.10 


± 


0.20 


0.02 


± 


0.07 


Geneva photometry 2 ' 3 


41* 


45 


± 


75 


0.00 


± 


0.35 


0.10 


± 


0.09 


Stromgren photometry 3 


• 4 50* 


190 


±110 


0.25 


± 


0.70 


-0.10 


± 


0.66 



1 Santos et al. (2000, 2001), Gonzalez et al. (1997-2001), Fuhrmann (1998), Edvardsson et al. 
(1993) 

2 Calibration by Kiinzli et al. (1997) 

3 Data from the GCPD (Mermilliod et al. 1997) 
4 Calibration by Napiwotzki et al. (1993) 

"The following solar type stars are included in addition to that listed in Tables 4 and 5: 
HD 16141, HD 18445, HD 22049, HD 29587, HD 95128, HD 110833, HD 112758, HD 140913, 
HD 168746, HD 178911, HD 202206. 

4.2. Fe abundances 

The parameters and Fe abundances determined spectroscopically in this work 
are given in Tables 4 and 5. Figure 10 shows the Fe abundances as a function of 
effective temperature for the two different groups of stars. For the calculation of 
the model atmospheres we used ODFs with solar abundances scaled according 
to the derived Fe abundances. We found this to be particularly important for 
overabundant stars, for which the abundances were underestimated when ODFs 
with solar abundances were used. 
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4. Parameters and Fe abundances for CGP dwarfs. 



HD 


T e s 


lo 


gff 


^micro 


[F 






8574 


6100 


4. 


30 




45 


-0. 


04±0. 


.05 


9826 


6200 


4. 


.40 


1 


.60 


0. 


13±0. 


.07 


10697 


5650 


4. 


.05 


n 


90 


0. 


16±0 


.06 


12661 


5700 


4. 


.10 


n 


50 


0. 


39±0. 


.07 


19994 


6150 


4. 


30 


l 


.50 


0. 


23±0. 


.07 


28185 


5700 


4. 


.20 


n 


80 


0. 


24±0 


,06 


33636 


6050 


4. 


.55 


l 


.25 


-0. 


,11±0. 


,05 


37124 


5700 


4, 


.50 


l 


.10 


-0 


,38±0. 


.05 


38529 


5750 


4. 


05 


l 


50 


0. 


48±0. 


12 


50554 


6050 


4. 


.50 


i 


,20 


-0. 


,04±0 


,04 


52265 


6100 


4. 


.40 


l 


.30 


0. 


15±0. 


.05 


68988 


6000 


4. 


.45 


l 


.35 


0. 


36±0. 


,06 


74156 


6050 


4, 


.35 


l 


.25 


0. 


.09±0. 


.05 


75732 


5500 


4. 


.40 


i 


00 


0. 


,55±0 


.16 


80606 


5700 


4. 


.40 


0. 


.90 


0. 


46±0. 


.07 


82943 


5900 


4. 


.40 


0. 


.75 


0. 


24±0 


.04 


89744 


6300 


4. 


.40 


i 


,80 


0. 


22±0. 


.08 


92788 


5700 


4, 


.20 


0. 


.50 


0. 


.30±0. 


.06 


106252 


5890 


4. 


.40 


1 


1 


-0 


,10±0 


,04 


114762 


6000 


4. 


.40 


1 


50 


-0 


,78±0. 


.06 


117176 


5600 


4. 


.10 


1 


.00 


-0 


,05±0. 


.03 


120136 


6600 


4. 


.70 


1 


.90 


0. 


37±0 


,12 


130322 


5400 


4. 


.40 


n 


00 


0. 


,10±0. 


.06 


134987 


5720 


4. 


25 


0. 


.70 


0. 


35±0. 


.11 


136118 


6250 


4. 


.45 


i 


60 

■ \J\J 


-0 


,03±0. 


.07 


141937 


5900 


4. 


.45 


0. 


.90 


0. 


12±0 


.04 


143761 


5900 


4. 


.40 


1 


.25 


-0 


.25±0. 


.04 


145675 


5500 


4, 


.30 


1 


.10 


0. 


.59±0. 


.15 


168443 


5600 


4. 


.10 


n 


80 


0. 


09±0 


.04 


169830 


6300 


4. 


.40 


i 


60 


0. 


09±0. 


.06 


177830 


5000 


3. 


.70 





.90 


0. 


61±0 


,14 


178911 


6050 


4. 


.50 


l 


.10 


0. 


15±0. 


,06 


179949 


6200 


4, 


.50 


l 


.20 


0. 


20±0 


.06 


186427 


5800 


4. 


.40 


0. 


.95 


0. 


06±0. 


.03 


187123 


5800 


4, 


.35 


0. 


.90 


0. 


08±0 


.04 


190228 


5360 


3. 


.90 





.90 


-0 


17±0. 


.04 


192263 


5100 


4. 


.45 


0. 


.00 


0. 


11±0. 


.10 


195019 


5830 


4. 


.30 


1 


.05 


0. 


05±0. 


.04 


209458 


6100 


4. 


.50 


1. 


.30 


-0 


,02±0 


,05 


210277 


5700 


4. 


.40 


1. 


.20 


0. 


.27±0. 


.06 


217014 


5750 


4. 


.25 


0. 


.70 


0. 


,24±0 


,05 


217107 


5750 


4. 


35 


1 


.15 


0. 


39±0 


,06 


222582 


5800 


4. 


.40 





.80 


0. 


03±0 


.04 


BD -10 3166 


5550 


4. 


.45 


1 


.00 


0. 


51±0. 


,10 
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Table 5. Parameters and Fe abundances for no-CGP dwarfs. 



rilJ 




lo, 


gg 


^micro 


rural 




166 


5550 


4. 


50 





.80 


0.13±0. 


05 


4628 


5150 


4. 


60 





,80 


-0.21±0. 


.08 


10476 


5200 


4. 


35 





.00 


0.03±0. 


.07 


12235 


6100 


4. 


40 


1 


.60 


0.35±0. 


.11 


16160 


5100 


4. 


.55 





.60 


-0.03±0. 


.14 


16895 


6500 


4. 


.70 


1. 


.70 


0.07±0. 


.08 


22484 


6050 


4, 


.30 


1 


,55 


-0.07±0. 


.09 


26965 


5300 


4. 


.55 





.70 


-0.26±0. 


.08 


32147 


5400 


4. 


.65 


1 


.95 


0.50±0. 


30 


48682 


6200 


4. 


.60 


1 


.30 


0.13±0. 


.13 


50281 


5100 


4. 


.50 


1 


.40 


0.00±0. 


13 


76151 


5700 


4. 


35 





.65 


0.10±0. 


05 


84737 


5950 


4. 


.30 


1 


,30 


0.07±0, 


.04 


126053 


5650 


4. 


.40 





.65 


-0.44±0. 


.04 


149661 


5300 


4. 


.40 





,00 


0.15±0. 


13 


157214 


5650 


4. 


35 


0. 


.60 


-0.44±0. 


.06 


166620 


5200 


4. 


.50 


0. 


,50 


-0.09±0. 


.10 


170657 


5200 


4. 


.55 


0. 


,60 


-0.11±0. 


.08 


185144 


5400 


4. 


.50 


1 


.00 


-0.20±0. 


.06 


186408 


5780 


4, 


.35 


0. 


.85 


0.09±0, 


.04 


201091 
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Figure 10. Fe- abundances (relative to the Sun) of CGP and no-CGP 
dwarfs. 
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5. Conclusions 

At this point of time we cannot draw any hard conclusions, because 

• the sample sizes are limited; 

• the analysis is still incomplete; 

• the samples could be affected by selection bias. 
However, Figure 10 indicates that 

• in contrast to the no-CGP sample, there seems to be a lack of cool stars 
{T e s < 5400 K) in the CGP sample. Does this mean a bias in the sample 
selection? 

• there is a "clump" of metal-rich CGP stars at solar temperature. Is this 
a lack of metal poor giant planet hosts or a bias - either from selection or 
from a lack of metal poor solar type stars (with or without planets)? 

• there seems to be no correlation between metallicity and planet hosting at 
hotter temperatures (T e g > 6000 K). 

Planned future work includes: 

• Abundance analysis of a second comparison sample: Stars which have 
spectral types like those of the dwarfs with (CG) planets selected randomly 
from the Bright Star Catalog (and Supplement); 

• Abundance analysis of "Very Strong Lined" dwarfs (Eggen 1978); 

• Abundance determination for elements other than Fe, for all stars of the 
four samples. 
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